#Exercice 1

setwd("C:/Documents and Settings/A276/Bureau")
thc <- read.table("thc.txt", sep=" ", col.names=c("Groupe","Activ"))

#1.
A <- thc$Groupe == "A"
ActivA <- thc$Activ[A]
moyA <- mean(ActivA)
varA <- var(ActivA)
B <- thc$Groupe == "B"
ActivB <- thc$Activ[B]
moyB <- mean(ActivB)
varB <- var(ActivB)
C <- thc$Groupe == "C"
ActivC <- thc$Activ[C]
moyC <- mean(ActivC)
varC <- var(ActivC)
D <- thc$Groupe == "D"
ActivD <- thc$Activ[D]
moyD <- mean(ActivD)
varD <- var(ActivD)
E <- thc$Groupe == "E"
ActivE <- thc$Activ[E]
moyE <- mean(ActivE)
varE <- var(ActivE)

hist(ActivA)
hist(ActivB)
hist(ActivC)
hist(ActivD)
hist(ActivE)

#3.
bartlett.test(list(ActivA,ActivB,ActivC,ActivD,ActivE))

#4.
moyActiv <- mean(thc$Activ)
dl1 <- 5-1
dl2 <- 47-5
totdl <- dl1+dl2
SSB <- (10*moyA^2) + (10*moyB^2) + (9*moyC^2) + (8*moyD^2) + (10*moyE^2) - 47*(moyActiv^2)
SSW <- ((10-1)*varA) + ((10-1)*varB) + ((9-1)*varC) + ((8-1)*varD) + ((10-1)*varE)
SST <- SSB + SSW
MSB <- SSB/dl1
MSW <- SSW/dl2
F <- MSB/MSW
DL <- c(dl1,dl2,totdl)
Somme <- c(SSB,SSW,0)
Moyenne <- c(MSB,MSW,0)
Rapport <- c(F,0,0)
data.frame(DL,Somme,Moyenne,Rapport)

analyse <- aov(Activ~Groupe,data=thc)
summary(analyse)
